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I. INTRODUCTION 

The interaction of a J/ip meson with strongly interacting matter is to date still a contro- 
versial subject. While the production of J /ip can be understood within perturbative QCD 
due to the large mass of the charm quark, its further interaction with surrounding matter 
is essentially soft in nature and as such not treatable perturbatively. The knowledge of 
hadronic interactions, as well as their modifications at finite temperature and density, is 
however necessary for a proper understanding of ultrarelativistic nucleus-nucleus collisions, 
especially in view of a possible transition from hadronic to quark matter |jl| . The suppression 
of J/ip was initially proposed as a signal for a quark-gluon plasma 0]. Such a suppression 
was observed by NA38 0]. However, the data can be described by a variety of models on a 
phenomenological basis, both in a plasma [Q and in a conventional hadronic scenario [§|,§|. 
Thus, the question of the significance of the J /ip signal remains as yet undecided. 

Plasma formation is not expected to occur in hadron-nucleus (hA) collisions. However, 
data taken in hA collisions already show a considerable reduction of J/^ production at low 
xp (the region where J /if) is also measured in nucleus-nucleus (AB) collisions) as compared 
to proton-proton {pp) [|^,|^ . The suppression pattern in both hA and AB collisions is found 
to be consistently described by a phenomenological absorption cross section of af^s ~ 5-7 
mb§. 

On the other hand, it was recently argued that, due to the smallness of the heavy quark- 
antiquark system (QQ), a gluon needs to be sufficiently hard in order to resolve this pair, 
Ql > 1/''^^ Ill0|- Oiily deconfined matter at temperatures beyond the phase transition 
temperature Tc contains sufficiently many hard gluons to cause a J/ip suppression of the 
observed magnitude. Matter in the form of hadrons does not provide enough hard gluons, 
and consequently a value of al^^^ ~ 5-7 mb has been regarded as unrealistic. This obvious 
contradiction is one example for the need of an understanding of hadron-hadron interactions 
on a more fundamental level. 

With the present work, we aim to provide a step towards filling this gap. The approach 
we use is the description of hadrons as bound states of quarks. It allows one to consistently 
account for substructure effects. A full treatment of the hadron-hadron interaction is to date 
not possible due to the non-perturbative character of QCD in this region. A microscopic 
calculation of the J/ijj breakup process by impact ionization has been performed within 
perturbative QCD in Ref. fllj for a dense partonic environment and in Ref. |jTO| for a hadronic 



medium. However, in these approaches non-perturbative correlations in the final state 
(charmed hadrons) have been neglected, i.e. only the breakup of a J/ip into free charm quarks 
was considered. As lattice gauge simulations of QCD suggest ||T^-[l4|, hadronic correlations 
persist even for temperatures well above the deconfinement transition. Therefore, effective 
approaches to J/ip dissociation in the non-perturbative domain of strongly correlated quark 
matter consider this process as a quark exchange (string-flip) process []T3|. The role of 



quark exchange processes in hadron-hadron interactions has been investigated in several 
approaches [|T^-^. Recently, a systematic analysis of quark exchange contributions to the 
meson-meson interaction has been given in Refs. |21|J2^ within a diagrammatic technique. 
These approaches use a non-relativistic quark potential model to describe mesons as bound 
states. They have been applied to the elastic scattering of light mesons. When translating 



the diagrams into the language of Green functions |22], a generalization to flnite temperatures 



and densities as well as to a relativistic effective meson theory is possible. 

In the present work, we extend this diagrammatic technique to the calculation of the 
cross section for the inelastic reaction {QQ) + {qq) -^ (Qq) + (qQ), where Q and q stand for 
heavy and light quarks, respectively. We consider the process of charmonium dissociation 
by inelastic collisions with light mesons. In particular, we calculate the cross section of the 
reaction J/ip + tt — > D* + Z) as a function of the relative kinetic energy of the mesons and 
address its application to the analysis of the kinetics of charmonium dissociation in heavy 
ion collisions. 

The section ^ gives the general formalism for inelastic meson-meson scattering. In 
section [TITI , the special case of charmonium dissociation is treated within this formalism, 



and cross sections for the main processes are calculated. These are then used in section |V 
to study the absorption in a pion gas. In section |V] we discuss the situation encountered in 
the experiment. 

II. QUARK EXCHANGE CONTRIBUTION TO MESON-MESON SCATTERING 

In this section, we present the formalism for the calculation of the cross section of quark 



exchange processes between mesons. We mainly follow the notation of Ref . |^ . We consider 
the two-meson scattering process A{aa) + B{bb) — >• C{ab) +D{ba), where the interchange 
of the quark content (in brackets) corresponds to a flavor rearrangement. This process 
dominates the cross section behavior at low relative energies of the mesons, while at higher 
energies the additional production of light qq pairs sets in which is not contained in the 
present approach. 

The differential cross section for the process i — > / is given by 

where P{s) is the relative three-momentum of incoming particles in their center of mass 
frame. For the relation between P and s, t see App. 0. The indices i and / stand for the 
initial and final two-meson states. The central problem is the calculation of the relativistic 
invariant matrix element Aifi{s,t). For this, effective theories have to be used in the low 
energy domain, where perturbative QCD is not applicable. One specific property of hadron- 
hadron scattering is the color neutrality of asymptotic states, such that a single one gluon 
exchange between hadrons is forbidden. The quark exchange process, however, is possible 



and the matrix element reads in Born approximation |21,22 



Mf, = M (^^^^ Hab,cd ^''^'') , (2) 



with the meson-meson interaction Hamiltonian Hab,cd and a product ansatz for the in- 
coming (outgoing) two-meson states formed by the mesons A,B {C,D). Four-momentum 
conservation is implemented in this matrix element. The normalization factor A/" is needed 
in order to get the correct form of A^/i from the nonrelativistic transition matrix element. 
With our convention, it reads 

A^=-^ n \/2E^^o, (3) 

"0 i=A...D 



where Qq is the norinahzation volume of the states \1/, which is set in the following to unity, 
and Ei = Jmf +pf. The calculation of M.fi is performed in the center of mass frame of 
the mesons A and B. The resulting differential cross section da/dt is expressed in terms 
of Mandelstam variables, that is, in Lorentz invariant form. The total cross section for 
scattering into channel / is obtained by integrating over t 

where t+(t_) is the maximal (minimal) possible momentum transfer t. The t integration 
can be transformed into an integration over z = cos'i9(P, P'), where '(9(P,P') is the angle 
between the relative momenta P and P' of incoming and outgoing mesons, respectively. For 
nonidentical particles the following relation holds: 






j'^dz\Mf.{P{s),P'{s),z)\'. (5) 



A. The quark exchange Hamiltonian 



It has been shown in Refs. ||21| , |22| , that in the quark potential model the Hamiltonian 
of the quark exchange process in meson-meson scattering can be represented as an effective 
two-quark interaction followed by a quark interchange between the mesons. The result of 
the diagrammatic analysis of all topological inequivalent contributions to the quark exchange 
matrix element is shown in Fig. |I|. According to Eq. (H) the matrix element of the quark 
exchange Hamiltonian in the four quark basis reads 

(^a,a,b,b\HAB,CD\c,c,d,dj = J2 [a,a,b,b\H-j\c,c,d,dJ. (6) 

For illustration we give the first of these four terms 

a,a,b,b\H^l\c,c,d,d) = ^^ (^a,b\H^\a\Vj {a,b\a\b') (a\a\b\V\c,c,d,d) , (7) 

a'.ii' 
6', 6' 

where a . . .d' denote three-momentum, spin, flavor and color quantum numbers of the quark 
or antiquark (a = {p^, s^, /a, Ca}). The last bracket selects those contributions in the sum 
over all quark quantum numbers, which match to the final state two-meson wave function. 
The other terms in Eq. @ are obtained in an analogous manner, where the interaction acts 
between the particles i and j. 

In the quark potential model, the two-quark interaction of a meson is given by the 
interaction Hamiltonian H^ of Fermi-Breit type. The same interaction is assumed to act also 
between the quarks of different mesons. It consists of the usual kinetic term, a nonrelativistic 
potential {H^) and relativistic corrections arising from spin-spin (H^^) and spin-orbital 
(H^^) interaction, a tensor interaction {H^) and a spin independent term (H^^), see ||2^ 
for a review. In the present work only S wave mesons are considered. Thus the two-quark 



interaction Haniiltonian H^ is a sum of the quark-quark potential H^ and the spin-spin 
interaction H^^ only, 



^,J\Hy,J')={^,J\H'\^\J') + {^,J\H''\^\J'). (8) 

For the reason of mathematical tractability of the matrix element (0) we choose an effective 
Gaussian ansatz for the orbital part of interaction, which in momentum space reads 

a, -h\H^\a\ V) = -l-o(87rx)^/^e-2^(P"-P'^)^5i;:;f'->5iy-'5p„^p_^,p, +p, , (9) 

with parameters Vq and x for the depth and range. A Kronecker symbol with superscript 
S,F or C is understood to act in spin, flavor or color spaces, respectively, e.g. ^jj/^^ = 
^s ,s-,Sf ,f-,Sc ,c,- This potential does not account for confinement. The use of a nonconfining 
potential follows Ref. ||2l| and is justified as long as the wave functions become small in the 
vicinity of the edge of the potential. In particular for the mesons containing heavy quarks, 
this condition holds. The spin-spin interaction is taken in the standard form p3[ 



^'^1^""^' ^ / = 9^^^^^^5"^S^'^p.+p.pHp;> (10) 



32TTas 

lilTlj 

where rrii and nij are the constituent quark masses. This can be rewritten in terms of H^ 
(Eq. ^ since the identity operator in momentum space can be understood as a limit of the 



Gaussian potential for a; — ;► and Vq -^ (Svrx) 



-3/2 



B. The meson wave functions 

We decompose the mesonic wave functions into orbital ($), spin {xs), fiavor and color 
iXFc) parts, 

(a, a|^^) = $p^(Pa,Pa)X5(Sa,Sa)XFc(/a,/a,Ca,Ca). (11) 

Instead of finding the exact eigenfunction of the two particle Schrodinger equation we use 
trial Gaussian wave functions and find the best approximation by using the Ritz variational 
principle. The orbital part of IS" state wave function is given by 



'4A,^^/^ 



n 



$pJPa, Ps) = (27r)^/^ [^^\ exp -2Xap\ 5p.4,(p.+p.), (12) 



where A stands for the quantum numbers and P^ = Pa + Ps for the total momentum 
of meson A. The relative momentum of the quark and antiquark in the meson is p^ = 
VAPa ~ (1 "~ VA)Pa, where tja = rria/ima + m-a). The constant A^ is related to the mean 
squared meson radius via {r"^) ^ = 6A^. 

We calculate the matrix element in the center of mass frame of mesons A and B, where 
Pa + Pb = Pc + Pd = because of total momentum conservation. Let us introduce the 
notation 



P = P^ = -Pb, 

P' = Pc = -Pd. (13) 

The generalization of the wave functions to excited states is straightforward. If one would 
also consider P waves, the spin-orbit and tensor terms of the interaction Hamiltonian had 
to be taken into account additionally. The parameters of potential and wave functions are 
fitted to the masses of the vr, p, J /ipi i^', D and D* mesons, see App. 0. 

C. The transition matrix element 

According to the diagrammatic analysis of the contributions to the matrix element (^, 
there are four contributions to be evaluated, see Fig. |l]. The first two diagrams correspond 



to the so called capture diagrams of Ref. |21], since the interacting quarks are captured in 



one meson in the final state. The others represent the transfer diagrams. The additional 
diagrams that arise, if identical quarks are present in the considered process, have the same 
amplitude and thus can be accounted for by a factor 2. 

Since H^ is a sum of an orbital and an spin-spin interaction, the transition matrix 
element can be written as a superposition 

Mf^ = E -^K + M!f. (14) 

We make use of the product ansatz for the wave functions Eq. (|TT|) and calculate the con- 
tributions from the terms of Eq. ([T^ ) to the matrix element M.fi- Each of these factorizes 
into an orbital (/o), spin {Is) and flavor-color (/fc) part. H^ acts on the orbital part of 
the wave functions, and H^^ on the spin part. For the matrix element M.^^ we obtain 

Ml{P, P', z) = ^{P, PVodP^ P'^ ^) ^h ^Fcp (15) 

with 

lldP.P'.z) = ($^<l>.i|i^J|$^.$5>.), (16) 



Ih, = (X^xf Ixk^), (17) 

FC,ij 



Ifc, = {xicXfc\x%X^c)- (18) 



The calculation of the orbital, spin and flavor-color factors is explained in App. 0. Here 
we only give the result: 



M:-,{P, P', z) = -Af{P, P') ll,-,4c,a-bKaiew [- [<yiP' + «2P'' + asP'Pz)\, (19) 

M^P, P', z) = -Af{P, P') ll,,I^c,abKabexp [- ((3,P^ + f3,P'' + P^P'Pz)] . (20) 

The constants ai, ■ ■ ■ , Ps, K^i and Kab are fixed by the parameters of the potential (Vq, x) 
and the wave function (Aj) and are explained in Eqs. (|U9D-( PT^ . It can be shown that the 
M-ab diagram is obtained from M.al by interchanging mesons C and D and replacing z by 
—z in (|C5|). Thus, one can express the matrix element A^^^ in terms of M^i by exchanging 



rjc ^^ 1 — f]D and A^ ^^ Xd and inserting the corresponding spin and flavor-color factors 
Ig and Ip(j. The same relations are valid between Ai^h and Ai^i ■ P' is fixed by P due to 
energy conservation. 

The corresponding matrix elements from the spin-spin interaction are 



ss. 






P,P',z) 



U{P,PV'o%{P,P' 



tSS tSS 
^S,ij ^FC,ij ' 



(21) 



with 



tSS 



(P, P\ 



($ 



A^B 



•I'g,^^ 



'), 



jSS 

tSS 
^FC,ij 



{XsXsWij \XsXs)^ 



(22) 
(23) 
(24) 



In this case, the Hamiltonian acts on the spin part of the wave functions. The spin factors Jj 

and I§^ are given in Table |. Because of the l/{ma'mi) dependence, the spin-spin interaction 
HSS 

Ref. 



dominates the matrix element when the interacton of light mesons is considered. In 
Jlp, the orbital interaction H^ has been disregarded in the calculation of the ti^tt^ 



scattering phase shifts, see also the following subsection. For our present application to the 
charmonium dissociation process, the contribution of H^ plays the dominant role in the 
transition matrix element ([T^) and it will be examined in more detail in section [ITl . 



D. Elastic 7r+7r+ scattering 



In this paragraph, we give the instructive limiting case of four equal quark masses rua = 
. . . = 1712 = "^g- That is, we consider the scattering of identical spinless IS mesons with 
masses m, described hj Xa = ■ ■ ■ = Xd = X. In this case the absolute values of the incoming 
and outgoing relative momenta P and P' are equal. We multiply by a factor of 2 in order to 
take into account the diagrams which arise in addition to those containing distinguishable 
particles. In this case, the matrix elements get a transparent form. The spin-spin term reads 



A<^^(P,z) 



32TTas 
drrig 

-|-exp 



i-2Q'\., 



-XP'{l + z) 



-4:X/3P^ 

XP\l 



exp 



)]}■ 



(25) 



It has been shown in |lT5| , pT| , that the low energy scattering phase shifts of tt+tt"'" scattering 
can be well described by this matrix element. The minimal relativistic approach to quark 
exchange processes in hadron-hadron scattering has also proven successful in the description 
of Kit and KN scattering p^. We expect that in processes, where quark creation and 
annihilation is negligible, the presented approach will be applicable. 



III. CHARMONIUM DISSOCIATION 



In this section, we apply the presented formalism to the specific case of the breakup of 
J/ip when scattering on hadrons. We calculate the absorption cross section from the quark 



exchange process, which is a function of the relative kinetic energy of the two scattering 
mesons, and discuss the imphcations for reahstic physical situations. To demonstrate the 
importance of correlations in initial and final states, we also discuss the breakup reaction 
of J /ifj into free quarks. Here, the result corresponds to previous perturbative calculations 

Charmonium absorption processes in hadronic matter have been considered in several 
works, e.g. P,|25|J2B[|. Basic processes for charmonium dissociation in hadronic matter are 



a) J/i) + -K^ D{IS) + D{IS) Am > 0.643 GeV, 

b) J/iP + p^ D{1S) + D{1S) Am > -0.13 GeV, 

c) J/tfj + N -^A^ + D{1S) Am > 0.258 GeV. (26) 

Generically, we denote D^,D^ or D'^ as D, and correspondingly D for the antiparticles. 
D{1S) represents either D or D*. Note that the reaction J/i/j + tt —>■ D + D, without excited 
final states, is forbidden by angular momentum conservation. The reaction thresholds for the 
possible processes are given by the mass differences Am. All these reactions are examples 
of inelastic quark exchange processes among hadrons. In the following, we work out our 
formalism considering process a) which describes charmonium absorption in a pion gas. 
Other processes including higher meson states such as Xc and ip' can also be considered, see 
the discussion in the conclusions. Process c) describes J/ip absorption on nucleons and can 
be treated on a similar basis. 

A. J/ip absorption by pion impact 

We apply the approach given in section H to calculate the energy dependent cross section 
of the process a) by specifying the initial mesonic states A = J/ip [QQ), B = tc {qq) and the 
final state , C = D{1S) (qQ), D = D{IS) (Qq), where Q is the heavy charm quark and q 
the light u or d quark. 

In order to work out the details in a transparent way, we use Gaussian wave functions 
and a Gaussian shape for the interaction potential which binds the quark-antiquark pairs 
into mesons. With the parameters of App. p, we obtain a satisfying description of the 
relevant meson spectrum, see Table |IT|. The choice of the Gaussian class of functions has 
the advantage that the calculation of the cross sections can be performed analytically, which 
makes the results more transparent. Eq. ([T9|) and (|20| ) are now used to calculate the cross 
section for the charmonium dissociation reaction J/tp + vr — >■ D{1S) + D{1S). Due to the 
large charm mass, the spin-spin interaction is negligible and we keep only the potential 
interaction. Then, for each final state channel, four matrix elements have to be computed. 
The integral over z in Eq. (0) can be performed analytically with the result 



X I Kl^ exp [-2aiP2 - 2^2^^'^ 
+ir2,exp[-2/3iP2_2/?2P'2 



sinh(2a3P'P)^ 
^ 2a3P'P 
sinh(2/53P'P)' 



1 + 



WiP'P 



Here, the first term arises from tlie two capture diagrams. Tliey liave the same spin and 
flavor-color factor and differ only in the sign of z, thus Ai^i{z) = Aiab{—z) and consequently 
K^l = Kab- In an analogous manner, we have ^Aab{z) = Mab{~^) ^^^ ^ab = K-^i for the 
transfer diagram that corresponds to the second term of the equation above. The last term 
contains an interference of both processes. We have defined /J := /J^^ = I^■^^^ = Jj^^, = /J-^ 

and I^c ■= ^FCab = ^Fc,ab = -^FC,ab = -^FCab- ^rom the parameters of the po- 
tential model (|B§, we obtain the values al=1.37GeV"^ a2=1.22G'e1/-2, a3=0.059 GeV'^ 



/?i=0.717GeV-^ /?2=0.507GeV-^ and /?3=0.368 GeV'^ Inserting this result in (D, with 
P{s) and P'{s) from Eq. ([X5|), we obtain the cross section ct/j for a specific final state /. 
For the total J/ip breakup cross section due to pion impact, we sum the possible final state 
combinations of low lying D mesons to get 

4 

o^abs(s) = Ea/.(s), (28) 

/=i 

where s is the center of mass energy of the J/ip and vr. The resulting J/ip absorption cross 
section, which is a function of the relative kinetic energy of J/ip and vr in the cm. system, 
^rei'^ = s — (m^ + m^)^, is the central result of this section. We show it as a function of 
E!^^'^ in Fig. |^. Here, the parameter values from Eq. (B4) are used and all low-threshold 



processes according to Table |I| are included except the lowest DD channel which is forbidden 
by angular momentum conservation. 

The behavior of the cross section is characterized by a threshold at Sq = {mc + moY and 
a strong enhancement near this threshold, as well as an exponential fall-off towards higher 
energies. We obtain a peak value of about 15 mb at -E^e?** — ^ GeV. The absorption cross 
section is approximately described by the fit formula 

a'^Xs) = a,-{l- ^)' exp [-a{s - s,)] e{s - s,). (29) 

The fit parameters for different possible final states ctq, a and si are given in Table |I|. 

As we mentioned in the beginning, we do not consider the inelastic production of ad- 
ditional light qq pairs, which sets in at a threshold of y^ + m^ + mq. Therefore, the 
exponential decrease in Eq. (^) is not considered to be realistic in view of the additional 
final state channels opened beyond this energy. 



B. Phenomenology of hadron hadron cross sections 

The large value of the absorption cross section we obtained within our calculation is at 
first sight a rather unexpected result. However, what was calculated is the cross section of a 
preformed, full-size J /tl) on a vr. This is in most situations not realistic. In real life, the QQ 
pair expands from a small object at the creation vertex to its full size [p7| , p8| . The initial 



8 



size can be estimated to be {f'^)QQ ~ l/(2^^c) ~ 0.06 fm, as supported by charmonium 
hadroproduction and photoproduction experiments. A quantum mechanical treatment of 
the expanding QQ state which simultaneously interacts with hadrons gives a time scale of 



this expansion of t'^J^ = 0.85 fm in the QQ rest frame |^| 



In the present hA experiments, the kinematics are such that asymptotic J/ip^s are ob- 
served only at high momenta in the final state. Then, r^^ has to be multiplied by a rather 
large 7 factor. In other words, the J/tp is only formed far outside the nucleus. This has 
two consequences. Firstly, the QQ interacts inside the nucleus still as a correlated, but 
considerably small state. We investigate this situation by describing the initial QQ state 
with a wave function narrower than the one of the J/ip, which is done by changing the wave 
function parameter Xqq accordingly. What we find is a decrease of the breakup cross section 
with decreasing QQ size. More precisely. 



This confirms the the phenomenological Povh-Hiifner relation |30| of hadron-hadron cross 
sections. Therefore, in realistic experimental situations, the cross section of 15 mb is lowered 
according to the kinematical circumstances. Secondly, possible differences in the final state 
interaction of J /ip and ip' , as expected already from their difference in size, are delayed, and 
thus become invisible because the difference in their asymptotic states appears only after 
they have left the target nucleus. This is supported by the experimental observation of an 
identical depletion of J/ip and ip' in heavy nuclei 0. 

C. J/ijj breakup without final state correlations 

As mentioned in the introduction, it was argued recently within a perturbative approach 
that a J /ip breakup reaction via gluon exchange requires a relatively hard gluon in order 
to resolve the small QQ state [^. The result we obtain for the cross section. Fig. |^, 



is completely different from the cross section obtained in such a perturbative approach. 
However, the quark exchange process we considered is also very different from a gluon 
exchange and essentially nonperturbative in nature. We note in this context that the present 
treatment can be traced back to older works of Gunion, Brodsky and Blankenbecler on 
composite models of hadrons |[T^ . They showed that even in certain short-range interactions 
constituent exchange dominates over gluon exchange processes. 

In order to illustrate this important point in the context of our approach, we calculate 
the cross section for the breakup reaction of J /ip and vr into four asymptotically free quarks 
within our effective model. Instead of the Gaussian wave functions (|T2D we define the final 
state as plane waves, in momentum space representation 

$p'(Pc,Pc) = '^P',{pe+Pc)^(pc-Pc)/2,PC' 
•^-P' (Pd, Vd) = '^-P',(pd + pj)'^(pd-pj)/2,pc • (31) 

The final four quark state is defined by the momenta pc? Pd and P'. The state C contains 
the charm quarks Q, Q and D the light quarks g, q. In diagram Fig. ^ all spin and color 
states are degenerate in the final state, and the sum over spin, flavor and color quantum 
numbers gives a factor 1 for Ig and Ipc- 



As before, the spin-spin interaction is small and we consider the potential contribution 
to A4fi only. For Ai^Jf^ we obtain 

X ;<(pc + ?) - $o^(Pc - f ); 

X $f(PD + f)-$f(PD-f)], (32) 

with <l>^ and $^ from Eq. (0), A = {QQ), B = (qq) and Q = P - P'. ^f{s,pc,-pD) is 
given according to Eq. ( |A8|) by 



^'{s, PC, Pd) = ^{^' - K - <)"} {'" - (4K + Pc) - 4K + pD)"} ■ (33) 

Inserting this into Eq. (^ we obtain the cross section into one definite momentum configura- 
tion crj7''(s, Pc, Pd)- The total cross section a^'^^'^{s) of the process J/ip + 'iT^Q+Q + q+q 
is given by integrating over all possible relative momenta pc and po 

free/^N [ dPC f dpDhccf, ~ ~ N fnA\ 

^ (*) = / 7KZ^ I TTTZ^^n (s^PcPd)- (34) 



(27r)3i (27r)3 ^* 



The integration is restricted by energy conservation to < p|, < {yi/2— (mg+p^)^'^} 
and < Pq < (^/i/2 — m^)^ — mg. Our result of the "perturbative" breakup cross section 
as a function of -E^^* is shown in Fig. |^. It does not exhibit a peak close to threshold, but 
starts smoothly and increases monotonically with energy. This is analogous to what has 
been calculated in QCD perturbation theory [|r^ . 

At low relative energies, the cross section of this process is small and it does not contribute 
to the J/ip disintegration. The comparison of the two cross sections, shown in Figs. |^ and 
^, demonstrates the importance of the correlation of the quarks in the final state. It has the 
consequence of a strong enhancement close to threshold where the relative momenta of the 
outgoing quarks are small and correlations between them are most pronounced. 

We emphasize at this point that the quark exchange reaction into correlated final state 
mesons does not proceed via intermediate free quark states. Therefore, the only energy 
barrier encountered in this process is the reaction threshold, i.e. the mass difference be- 
tween initial and final state mesons. It is understood in our approach as the difference of 
the respective binding energies, which is overcome by kinetic energy of the initial mesons. 
However, we stress that no intermediate energy barrier is present in this nonperturbative 
approach. This has to be seen in contrast to a perturbative calculation, where such a bar- 
rier occurs and where a nonperturbative mechanism, such as a tunneling process, has to be 
invoked additionally. 



D. Inelastic cross sections in the strange sector 

We briefly look at the related processes involving strangeness instead of charm. The 
meson-meson reaction in this case is(f) + Tr^K + K, which is, however, experimentally not 
accessible. Instead, we look at the baryonic reactions in the strange sector corresponding to 
the ones relevant for charmonium. These are 
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a) K~ +p-^ K + X, 

b) K+ +p^K + X, (35) 

and a review of the data is given in [0,^ . The cross section for reaction a) exhibits a strong 
peak at threshold and a subsequent decrease with increasing energy, while the cross section 
for reaction b) increases monotonically from threshold. At energies far above threshold both 
cross sections reach the same asymptotic value. The data qualitatively show exactly the 
behavior we expect. Process a) is dominated by a simple quark exchange process as we 
considered before, for which we calculated a strong peak at threshold, while reaction b) 
requires a hard ss production process since K'^ contains an s quark, while an s is needed for 
the A. Therefore, reaction b) does not show an enhancement at threshold. Although only 
being qualitatively, this strongly supports the approach presented here. 

IV. DISSOCIATION KINETICS IN A PION GAS 

In this section, we consider the relaxation of the charmonium fluctuation by string-flip 
processes in a dense hadronic medium such as the pion gas produced in a high energy 
nucleus-nucleus collision. To obtain the suppression of the bound QQ states, we fold the 
energy dependent absorption cross section calculated in the previous section with a thermal 
pion distribution which is chosen in a way to describe the pion multiplicity and shape of the 
rapidity dependence in the same reactions where the J /ip is measured as well. 



The time evolution of the J/ijj distribution is described by the Boltzmann equation p3|j3^ 



— /^(r, p^, t) + -^Vf^ir, p^, t) = -/^(r, p^, t) j 7^^/^(P^' "^^ ^) 

XO-abs[s(PV,P7r)]j(PV"P'r)5 (36) 

where s(p^, p^) is the center of mass energy and 



m'^ml 



'<"*■ P-' = e,(p,)eAp,) ''" 



is the flux of pions in the rest frame of the J/ip (see App. ^). Due to the small number 
oi QQ pairs, the inverse process of J/ip production in DD scattering is neglected, and the 
influence of the considered reaction on the pion distribution is negligible. The solution of 
Eq. ( p6D for an initial J/ip distribution /^(r, p^,to) reads 



/^(r, Pv,, t) = /v.(r - v^t, p^, to) exp 



dt' 



d^Pn 



to J (27r)3 



x/^(r - V^(t - t'), p^, tVabs[s(Pv>, P7r)]j(p,/., Ptt) 



(38) 



We are interested in the time evolution of the total number of J/ip^s resulting from 



the absorption by the breakup process considered in section HI. For a qualitative estimate 



we consider the survival probability of a J/ijj in a uniform thermal pion gas. In this case 
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the integration over the space coordinate r can be performed and the resulting momentum 
distribution of meson i is given by 

mipi^t) = J dhMr,pi,t). (39) 

The Boltzmann equation (|36|) simphfies to 



where the relaxation time t(p^, t) is defined in the rest frame of the pion gas by 

T{p^,ty^ = (crabs^^re/)„^P7r(i), (41) 

with the pion density 

pAt)=l^n^iPn,t). (42) 

The brackets denote the average over the pion distribution n^r which may, in general, be 
time dependent. 

1 f dp 

Xaabs[s(pv>,P7r)]j(PV',Pvr)- (43) 

In order to give a quantitative estimate of the relaxation time for the J/tp distribution 
in a dense pion gas, we consider the specific example of a thermal pion distribution in 
equilibrium as given by the Bose distribution 

U{E^, T) = 3 (exp [{E^ - fi)/T] - 1)-' , (44) 

where the factor 3 stands for the pion multiplicity. The temperature T and chemical potential 
fi may be chosen as time dependent for modeling the evolution of density and energy density 
in nucleus-nucleus collisions. For the chemical potential of pions we use the value /i^ = 
126 MeV with which the experimental heavy ion data can be well reproduced |3^. The 



temperature range from 120 to 210 MeV corresponds to pion densities from 0.22 to 0.84 
fm"^. 

We show our result for the thermal averaged cross section (crabs^rei)^ ^^^ different tem- 
peratures and momenta of the J/ip relative to the pion gas center of mass in Fig. ^. This 
quantity gives the mean capability of a pion to dissolve a J/ip which is moving with mo- 
mentum p^ through the pion gas. For low momenta of the J/ijJ, it corresponds to a small 
cross section since the relative energy exceeds the reaction threshold only in few collisions. 
The cross section then rises with increasing momentum of the J/ip. For all values of p^, the 
absorption cross section increases with increasing temperature of the pion gas. The resulting 
mean life time of a J/ip moving through a pion gas, defined in Eq.(36), is plotted in Fig. |^ 
as a function of the pion temperature for different J ftp momenta. In comparison to previous 
assumptions or phenomenological calculations |P^-|5B| , we find a rather strong absorption of 
the J/ip. This is caused by the enhancement of the energy dependent cross section when 
quark- antiquark correlations in the final state are taken into account. 
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V. DISCUSSION OF RESULTS IN VIEW OF THE EXPERIMENTS 

A. J /'ijj absorption in a pion gas 

For a qualitative discussion of this result, we compare the calculated relaxation time r 
with the mean life time of a hadronic fireball which is measured e.g. by interferometry in 
the NA35 experiment and found to be in the range of 5 fm/c for freeze-out temperatures 
T < ISOMeV 1^. This means that the relaxation time is of the same size as the lifetime 



of the fireball. From our results, we conclude that J/ip dissociation in a dense pion gas is 
capable of producing a rather large absorption. In particular, it is large enough to describe 
the J/ip suppression as observed by NA38. 

However, so far we have discussed an idealistic situation which is not met in heavy ion 
collisions. First, the assumption of an equilibrium pion gas without baryons and resonances 
is not realistic. Second, in the early stage of the collision, the densities are so high that 
a description by a free pion gas is not appropriate, and it is uncertain at which time such 
a description becomes valid. Third, as we discussed in detail in the previous section, the 
J/ip cannot be regarded as a fully developed object from the very beginning. For the latter 
two reasons, we overestimate the contribution of collisions with pions to the J /ip absorption. 
It is therefore interesting to analyze the variation of the effective cross section (cTabs'^rei) 
again as a function of the radius of the QQ wave function {t'^)qq- For the free reaction 
QQ + TT ^ £) + Z), we had obtained a cross section proportional to the radius squared of 
the QQ pair. We find that this relation still holds after averaging over the medium. The 
quantity (cTabs'^re/) is shown in Fig. ^ for different temperatures of the surrounding pion gas 
as a function of {r^)QQ- It is approximately proportional to the mean squared QQ radius. 
Due to the symmetry of the quark exchange process we obtain the relation 

,2\ /2\ 



{(T^hsVrel) OC (r ) (^ ) _ , (45) 



which can be regarded as a temperature- averaged version of the Povh-Hiifner relation [30 



It has been emphasized on the basis of hadron-nucleus experiments that there must be a 



considerable contribution to the absorption of J/ijj on nucleons [4C]. Within our formalism, 
absorption of J/tp on nucleons such as the reaction (22 c) can be treated in an analog 
manner. We expect qualitatively a similar result as we obtained for the breakup of J/^/^'s 
on pions. With a choice of the time-dependent size of the QQ system appropriate for the 
respective kinematics, the hadronic absorption of charmonium on nucleons should be well 
described by this approach, except for the very high momentum region {xp ~ 1) 0]. In 
nucleus-nucleus collisions, we have the additional absorption on the pion gas formed in the 
collision, on which we concentrated in this work. With the proper kinematics, it is reduced 
from the idealistic estimate above, but still to a value comparable to the absorption on 
nucleons. Therefore, our results show that the two hadronic absorption processes on pions 
and on nucleons, when taken together, are able to account for the observed J/tp suppression. 

B. Heavy ion beams and inverse kinematics 

In the near future, experiments will take data with the lead beam at CERN. It has 
been proposed to also study reactions in inverse kinematics, i.e. a heavy projectile on a light 
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target. In this setup, the J/ip formed in the hard colhsion will subsequently be taken over by 
the heavy projectile nucleus at a low relative momentum. Thus, it is an ideal tool to study 
the interaction of (almost) fully formed J/ip^s on nucleons. From perturbative calculations, 
no absorption is expected in this case due to the low density of hard gluons in nuclear matter. 
Our prediction differs from this. From our cross section, we expect a strong absorption of 
the J/ip as compared to the case of pA collisions. The limitations of the model we presented 
only allow us to make qualitative predictions. 

We consider the situation of the NA38 experiments with a beam energy of 200 GeV 
per nucleon, where dimuon pairs from J/ijj^s can be detected in the rapidity window of 
2.8 < y^ < 4.0. So we have in the mean y = 3.4. We compare the case of proton beam 
on lead target (a: pPb -^ "ipX) with a lead beam on a light target (b: Pbp -^ 'ipX). Both 
cases differ in the rapidity difference of the detected J/ip^s and the nucleons in the heavy 
ion. In the proton beam situation (a) we have a rapidity difference of Ay = 3.4, which 
corresponds to a Lorentz factor of 7^ = 15 and a c.m.s. energy of a nucleon and a J/ip 
of ^/s^ = 9.88 GeV. In the lead beam experiment with inverse kinematics (b), we have 
yb = 2.6, 7f, = 6.8 and ^/sb = 7.06 GeV. The A dependence of the J/ip production cross 
section for a time dependent absorption has been studied, among others, in Ref. |^ . 

We adopt a linear time evolution of the QQ radius before r^xp- From our result, Eq. 
(^) , follows an evolution of Uabs in terms of the Q Q proper time r as 

, , fffo(s)^ ioir<Texp 
^abs(r,s) = <^ ^-- (46) 

[ cro(s) for r > T^xp. 

Let us introduce the ratio of the J/tp cross sections in direct (a) and inverse (b) kine- 
matical processes as 

R = ^^^^^. (47) 

(^Ap^ip 



Then, following Ref. ^^ , we find for the quadratic time dependence of the absorption cross 
section a ratio of 



R = exp 



A fcro{Sa) cro{sb) 



47rr2_ V 72t;2 ^2^2 



exp 



(48) 



where w^ ~ f & ~ 1 are the relative velocities of the J/ip with respect to the scatterers. We 
assume that initial state modifications and color octet absorption [Q in both pA and Ap 
collisions occur at short time scales inside the nuclei such that these effects cancel when 
taking this ratio. Inserting the values given above and using an energy independent cross 
section of o"o(sa) = cro(sfe) = 5 mb, we obtain R = 1.22. Concluding, we expect a strong 
suppression of J/ip production in the inverse kinematical regime as compared to the lead 
target. 

An enhancement of the J/ip absorption cross section near threshold, as we obtained for 
absorption on pions, leads to an even stronger suppression in this case. Taking, for example, 
c"o('Sfe) = l-5cro(sa) gives R = 1.38. On the other hand, perturbative estimates suggest for 
this ratio a value of i? = 1. Therefore, the experiment carried out in inverse kinematics is 
clearly able to discriminate between the different models, and thus to indicate the dominant 
physical processes. 
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VI. CONCLUSIONS 

The aim of the present work was to estabhsh a formahsm in which the absorption of 
J/ip mesons on hadronic matter is treated in a microscopic approach. We describe mesons 
in a potential model, and consider quark exchange reactions between two mesons as the 
model for inelastic reactions such as open charm formation in a J/^/'-hadron collision. For 
an exploratory calculation, we have used a Fermi-Breit Hamiltonian and a Gaussian ansatz 
for the quark-antiquark wave functions. In this model the energies of bound states are 
fairly well reproduced. In addition, the model can still be treated analytically and gives 
transparent results which we have discussed in detail. 

In this work, we concentrated on the reaction J /ip + 7t —>■ D -\-D and calculated the cross 
section for the breakup reaction as a function of the relative energy of the colliding mesons. 
For comparison, we also considered the breakup reaction in free quark states. The result of 
the latter calculation is similar to what is obtained in the corresponding calculation in the 
framework of short distance QCD. The comparison of both our results, as shown in Figures 
^ and §, demonstrates the importance of correlations of the quarks in the final state. 

When considering correlated quarks, i.e. mesons, in the final state, we find an enhanced 
cross section at the reaction threshold, that is, for low relative momenta of the quarks in 
the final state. In general we are able to describe a correlated QQ pair before it propagates 
to become a fully developed J/ip- For both kinds of final states, we find an increase of the 
absorption cross section proportional to the square of the (time dependent) system size of 
the QQ state. 

With a view to the application of this result to heavy ion collisions, we have estimated 
the suppression of J/^'s in a dense pion gas in thermal equilibrium. Since this is an idealized 
situation, the results can only be considered qualitatively. They suggest that in the temper- 
ature region accessible to present experiments, non-perturbative (string-type) correlations 
in the final state (asymptotically, these are charmed hadrons) as described in the present 
approach are crucial for the understanding of the J/ip suppression pattern. The enhanced 
cross section that we obtain at the DD threshold is sufficient to explain the J/ip suppression 
in the NA38 experiment as due to absorption by pions and nucleons. This is in contrast 
to previous claims that J/ip suppression could only be related to the quark-gluon plasma 
formation fTH]. We pointed out that a measurement comparing pA with Ap, in inverse 



kinematics, could clearly distinguish between both physical pictures, since the one based on 
perturbative calculations predicts basically no suppression, while in the one presented here, 
a strong suppression is expected. 

The formalism that we have presented in this work is rather powerful. A straightforward 
extension can simply accommodate higher states, such as the charmonium states {xc, V'O i 
mesonic resonances {p,uj, . . .), and baryons (A^, A). We mentioned that similar results are 
expected for charmonium absorption on nucleons. Since quite a sizable contribution to the 
J/ip yield stems from x decays, also the x ~ '^ cross section should be calculated. 

An analysis of the ip' absorption cross section would be of particular interest, since the 
ratio ip' to J/ip yields in pA and AB collisions is supposed to be free of initial state effects. 
This ratio was measured recently by NA38 [Q, and it seems to be a much clearer probe of 



the state of matter than the J/^l> to continuum signal. 

As a further outlook, we mention also the possible application of the present approach 
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to the interaction of J/ip^s with deconfined quark matter, see e.g. |T^. For this apphcation 
the present nonrelativistic calculation of the matrix element A^/j has to be improved by 
evaluations within a relativistic potential model |^5|], where the effects of chiral symmetry 
restoration and quark deconfinement at finite temperature can be included. 
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APPENDIX A: RELATIVISTIC KINEMATICS 



We consider the two particle process A + B ^ C + D. The Mandelstam variables s, the 
center of mass energy of A and B, and t are given by 



s = {Va + VbY = {Ea + EbY - (Pa + Pb?, 
t={VA- Vcf = {Ea - Ecf - (Pa - Pcf- 

The flux j of particle A considered in the rest frame of B is 

VaVb 



(Al) 
(A2) 



J = Vrel 



EaE 



A^B 



with 



>-g+)(g-g_^ 

2 



(A3) 



s± = [mA ± mB 



Vrel is the velocity of A in the rest frame of J/ip 

A{mAmBY 



-"rel 



1- 



[s — m. 



m 



2 \2' 

b) 



(A4) 



By inserting Eq. ( |A1| ) into ( |A3| ) one obtains the form of Eq. (|37| ) for j, which depends on the 
three-momenta of both particles. The following relations between three-momenta (defined 
by Eq.([l3|)) and Mandelstam variables are valid in the center of mass frame of particles A 
and B: 



P\s) = 


4^ 


P'\s) = 


1 
4^ 



{ml + ml) 



{ml + ml) 



— Amj\m.l 



— Am^Wjl 



2P'P = 2P'Pcos^ 



t 



m. 



2 , {s + m\- ml){s + ml 
m^j H 



r/i 



D) 



2s 



(A5) 
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Ea{s)Eb{s) = — (s^ - {m\ - ml) j , 

Ec{s)En{s) = ^ («' - {ml - ml)') . (A6) 

If four identical mesons with mass m are considered, these formulae simplify 

P2 = p'2^i(s-4m2), 

2P'P = t~2m' + -. (A7) 

The kinematical factor A/'(P, P') of Eq. (0) is given in the center of mass frame by 



Af{s) = A^EaEbEcEd 
1 



-y (-5^ - i.m\ - m|)2) (s2 - {ml; - m%Y). (A8) 



APPENDIX B: FIT OF THE MESON SPECTRUM 

For an exact solution the mesonic wave functions have to be calculated by solving the 
Schrodinger equation in the rest frame (Pa = 0). We use Gaussian wave functions of the 
form of Eq. (O) as approximating test functions. Then $(^\ depends only on one parameter 



A, and is denoted by $f ). The best fit is found by using the Ritz' variational principle. 
The Schrodinger equation reads 

^|<(Pa,Pa)) = mAK(Pa,Pa)), 

i7 = m, + ms+p^/2m„ + p?/2ms + //^ (Bl) 

where m^ is the mass of meson A and H^ = H^ + H^^ is the interaction Hamiltonian given 
in Eq. (|^). Only H^^ gives a spin dependent contribution to the meson mass. So we get for 
the mass of the 15* state 

mA = m^ + m, + ^±^^L±^ -vJl + ^\ '" - (^ - 5a)AM^^ (B2) 

I'oXArfiam-a \ 2x ) 4 

where Sa is the total spin of the meson A. Under the assumption that the orbital component 
of the meson wave function is equal for different spin states we can eliminate the spin 
dependent term. Averaging over the spin and isospin states, the spin-spin contribution to 
the Hamiltonian cancels and we get an averaged mass m°^ for IS* states |^E'(A)) which depends 
on the wave function parameter A. For the ground state one has to fulfill approximately the 
conditions 

m^"(A) = {'^^{\)\H - H^^\^^{X)) 

= -mT\\) + \mi-\\) 

3(ma + mia] , , / A \ .^ x 

= ma + ma + ^ ^_^ 1+ , (B3a) 

lbXm,ama \ Ix j 
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and 

d 



dX, 



(^^(A)|/7-i7*^|^^(A)) 



= (B3b) 



for each of the quark pairings qq, qQ,Qq and QQ with appropriate parameters for quark 
masses rrii and wave function parameters Ajj. We determine the parameters in (D and (|12|) 
from the meson masses oi Ti, g,D,D*,ric and J/ip to 



mq = 1.84 GeV, Xqq = 0.755 GeV'^ 

nig = 0.34 GeV, Agg = 3.05 GeV-^ , . 

x = 1.47GeV-^ AQg = 2.1GeV-^ 
\/ = 1.24 GeV. 



Table |IIT| shows a comparison of the known meson masses with the calculated masses with 



the model parameters (B^). Since some masses are uncertain, the averages are only ap- 



proximately. The calculated root mean squared radii of the states are given in the last 
column. 



APPENDIX C: CALCULATION OF TRANSITION MATRIX ELEMENTS 

The Born matrix element Aifi has to be calculated by integrating over all internal quark 
variables. Here we demonstrate this for processes with exchange of antiquarks of mesons A 
and B. If identical quarks are involved, one has to add to each diagram the corresponding 
ones decorated with a fermion commutation operator (with negative sign, see |^). The 
interaction can be written as sum of individual interactions between the quarks a, a, b and 
b, each of them consisting of potential (H^) and spin-spin [H^'^) contribution. We only 
show the calculation of the diagram A4^i of Fig. |1|. Calculation of the other diagrams is 
analogous. 

1. Orbital factor of the matrix element 



First we consider the orbital factor of Eq. (|T^), 

Io,aliP^P')= (<^P<^Sl^r5l<^P'<^-P')- (CI) 

In the center of mass frame the mesons A, B, C and D are moving with three-momenta 
P, — P, P' and — P'. From momentum conservation follows Pa + Pl = Pc + Pc, Pa = Pd) ^^^ 
Pb = Pd- Using these relations one can substitute six of the eight quark variables and has 
to sum up over the two remaining p^ and Pc- 



^^,.^(P,P') 



E ('^p|PaPs)('^S|P;>P6)(PaP6|^^|PcPc) 



Pa.Pa. 
Pd'Pd 



X (Ps Pb I Pd Pc) (Pc Pc I $P' ) (Pd Pd I $ P' ) 
Y: $*^(Pa - {l-VA)PW{Pa -P- VBP)H''{pa " Pc) 

Pa,Pc 

X$^(p, - (l-77c)P')$''(Pa - P - r^i^P'). (C2) 



The sums over p^ and Pc are replaced by integrals according to 

d^Pi 



E 

iPi> 



i2ny 



(C3) 



In our model we use Gaussian wave functions and a Gaussian shape of the potential. In 
this case an analytical expression for the orbital overlap matrix element in the (IS*) + IS* — »• 
(IS*) + (IS*) process can be obtained, 

-2 {Aa(p„ - (1 - VA)Pf + \B{Va - P' - VsPf + X{pa ' Pc)'} 

-2 {XciVc - (1 - Vc)P'f + MVa - P - VdP'?} 
= -K,i exp [- (aiP2 + aaP'" + a^PP'z) . 

For the general case we give only the constant ai, 

2 {A^Aslr/A + r^B - 1)' + AaA^(1 - r]Af + ^A^DiTlA? 

+XbXcvI + AbAz5(1 - Vb? + X'cXd} (Xa + Xb + X'^ + Xnr\ (C6) 



X exp 
X exp 



(C4) 
(C5) 



ai 



with 1/X'c = 1/Ac + ^/x. The other constants have similar forms. 

In the special case of the quark rearrangement reaction (p6|3) in section |T| we have due 
to flavor conservation 



rrin 



rua = rric 



m- 



^Q, 



rUb = mi = m-c = rud = niq 



We introduce the notation rj 



mq 



mq+mq 

Ac = Ad 



The outgoing D mesons have the same radii, 



Then the parameters of the capture diagram /q ^^ simplify to 



(C7) 



(C8) 



"^ = a:U 



+ Agg) [X'q^ + Aggj + Ag^A^ 



Qq\ ' 



«2 = 7- {(Aq9 + AqJ (XqQff + Agg(l - Vlf) + XqqXqq] 
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a; 



Aq = Xqq + Aqg + Xgg + Agq, 

In a analogous way we obtain Ipab^ 



Qq ' 



ab- 



(C9) 
(CIO) 

(Cll) 



Kab exp [- (AP2 + /32P'' + I3^PP'-^ 



(C12) 



where in the considered process 



2 1 1 

^[4^00(^9 + ^Qg + 2a;) + ^Agg(AQQ + Aq, + 2x) + Aq^x] 



/32 



A^ 



V^XQQ{\q + Agg + 2x) + (1 - rif\q{XQQ + Aq^ + 2x) + AggAggX/AQq 



/^S = -r- [^AQQ(Ag, + Xqq) + (1 - r7)Ag5(AQQ + Aq,)] 

A/3 






(Aqq + Xqq + x)(Agg + Xqq + x)-x^ 
V.(^-^y{XqqXqqfl^=K,,. 



I^Q^ 



qi 



(C13) 
(CM) 

(C15) 



P' is determined by P from energy conservation. The matrix element can be rewritten in 
kinematical variables s and t, see App. ^ The orbital factors for the remaining diagrams 
loab ^^d Iq-^i are obtained by replacing z by —z in (|C5| ) and (|C12| ) because the orbital 
wave functions of mesons C and D are identical. 



2. Spin factor 

The spin wave functions for spin singlets and spin triplets are given by Clebsch-Gordan 
coefficients 



Sa,S 



Xs )— Z^\^a-iSg^,Sa-iS^\SA-iSj^\Sa-iS^,Sa-iS-^). 






II = {XsXs\XsXs) 

= E xi(s.,s,)x|(s.s,)5i^5^^i5^fj5fMxS(SeS.)x!,(s.s-,) 



(C16) 
(C17) 



a,a,b 
...,d,d 



The different values of this factor are summarized in Table |. The spins of initial mesons 
and the sum of them are in the head of the table and of the reaction products in the first 
column. The spin factor in the potential interaction term of the transition matrix element 
between these states can be read from this table. For the If "^ term one has to multiply these 
factors by the numbers given in the right two columns for the diagrams Aial ^'^d Aiab- The 



factors for Aiab and Ai^l diagrams are obtained in a similar way, see [^ 
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3. Flavor-color factor 



The color singlet wave function for mesons is 

X^c = ^j: C (C18) 

and the flavor and color give only a conibinatoric factor in the Born matrix element. From 
the color component we get an overall factor of 1/3 for M-al ^-^d J^ab and -1/3 for M.ab 
and M-ai. This accounts for the fact that only a third of all quark pairs, those of identical 
color, are able to produce a color singlet in the quark exchange process. This result differs 
by a factor 4/3 from Ref. |2^. This factor is included in the effective form of our model 
interaction. The light quark flavors u, d are assumed to be degenerate. 



'FC 



{xicXFc\x%XFc) 

xUara)xfcibmTS^ffSi:;r6\,r'xUc-c)x'i^cidd). (C19) 



,B (hh\K^F,C)^(F,C)<^(F,C)^(F,C)^^A l^^\^,C 

a,a,b 
...,d,d 



From the Kronecker deltas follow that only the quark line diagrams shown in Fig. |^ give 
nonzero contributions to the matrix elements. All other possible diagrams are forbidden in 
our case of four different quarks. I^'" gives ±1/3 if flavor conservation is fulfilled and zero 
if not. 
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TABLES 

TABLE L 
Spin factors /J of the transition matrix elements for different spin states of the initial {A, B) 
and final (C, D) mesons. For the I§^ factors one has to multiply each line by the factors in 
the respective column /|^j or /ff^. 

TABLE IL 
The spin factors for different final channels / = {C, D) in the quark exchange process J /ip + 
TT -^ C + D. They follow directly from Table |. ctq, sq and a are parameters for the fit formula 
Eq. (HD. 



TABLE in. 
Meson mass spectrum according to the formula Eq. ( p3|) with the fitted model parameters 
( [B4| ) in comparison with the spin averaged experimental masses. In the last column, the 
root mean squared radii of the mesonic state are given. 
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FIGURES 

FIG. 1. 
Contributions to the quark exchange matrix element Ai/i in the Born approximation: 
M-ab^-^'ab (capture) and A^ab, -^afe (transfer). Each interaction line represents the sum 
of potential and spin-spin interaction. 

FIG. 2. 

Cross section for different channels of inelastic rearrangement reactions of J /ip and vr into D 
mesons. 

FIG. 3. 
Diagram for the disintegration reaction of J/^ and tt into four free quarks. 

FIG. 4. 
Cross section for the reaction J /ip + 'K^Q+Q + q + q. This result is comparable to the 
perturbative calculation of Kharzeev and Satz [0 . 

FIG. 5. 
Thermal averaged cross section {aVrei)T ^^ J 1'^ + ^ ~^ ^i^^) + D{1S) for a J/ip moving 
with momentum |p^| through a gas of pions with a chemical potential fi = 126 MeV at 
different temperatures T. 

FIG. 6. 
Mean life time r of a J/ip in a pion gas as a function of the temperature T for different 
momenta |p^| of the J/ip with respect to the pion gas center of mass. 

FIG. 7. 
Dependence of the thermal averaged cross section {aVrei)T o^i the mean squared radius of the 
QQ wave function for a, QQ pair at rest in an equilibrium pion gas with chemical potential 
/i = 126 MeV at different temperatures T. The nearly linear dependence for small ygg) 
corresponds to the Povh-Hiifner relation [RO]. 
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TABLE I: Martins et al., Quark exchange 



quarks 


2^+%i, 


meson 


m[GeV] 


me%[Gey] 


m'"'[Gey] 


{r-Y" 


qq 


^15o 


TT 
P 


0.136 
0.771 


0.612 


0.61 


0.84 


QQ 


^15o 




2.980 
3.097 


3.068 


3.07 


0.42 


QQ 


HP2 


XcO 
Xcl 
Xc2 


3.415 
3.511 
3.556 


3.525 


3.43 


0.54 


QQ 


^25i 


i'{2S) 


3.686 


W3.64 


3.69 


0.64 


Qq,qQ 


'l5o 


D 


1.868 
2.009 


1.974 


1.94 


0.70 


Qq,qQ 


HP^ 
'IP2 




2.424 
2.459 


W2.43 


2.38 


0.90 



TABLE II: Martins et al., Quark exchange 



/ 


{C,D) 


i/jr 


(To [mb] 


so[GeV"] 


si[GeV2] 


a[GeV-^] 


1 


(D,D) 








14.0 


- 


- 


2 


(D*,D) 


1/4 


4.33 


15.05 


24.43 


0.549 


3 


{D,D*) 


1/4 


4.33 


15.05 


24.43 


0.549 


4 


{D*,D*) 


1/2 


7.42 


16.2 


25.2 


0.621 



TABLE III: Martins et al., Quark exchange 
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FIG. 1: Martins et al., Quark exchange 
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FIG. 2: Martins et al., Quark exchange 
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FIG. 3: Martins et al., Quark exchange 
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FIG. 4: Martins et al., Quark exchange 
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FIG. 5: Martins et al., Quark exchange 






11 =126 Me V 




125 0.15 0.175 0.2 0.225 0.25 0.275 

T (GeV) 



FIG. 6: Martins et al., Quark exchange 
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FIG. 7: Martins et al., Quark exchange 



